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Abstract 

Unlike ground-based interferometric gravitational wave detectors, large space-based systems will 
not be rigid structures. When the end-stations of the laser interferometer are freely flying space- 
craft, the armlengths will change due to variations in the spacecraft positions along their orbital 

(N 

trajectories, so the precise equality of the arms that is required in a laboratory interferometer 

O ' 

to cancel laser phase noise is not possible. However, using a method discovered by Tinto and 
(3JT), Armstrong, a signal can be constructed in which laser phase noise exactly cancels out, even in an 

unequal arm interferometer. We examine the case where the ratio of the armlengths is a variable 
parameter, and compute the averaged gravitational wave transfer function as a function of that pa- 
rameter. Example sensitivity curve calculations are presented for the expected design parameters 
of the proposed LISA interferometer, comparing it to a similar instrument with one arm shortened 
by a factor of 100, showing how the ratio of the armlengths will affect the overall sensitivity of the 
instrument. 
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I. INTRODUCTION 



One of the differences between laboratory and space laser interferometer gravitational 
wave detectors is that, in the laboratory, the two arms of the interferometer that is used 
to detect changes in the spacetime geometry are maintained at precisely equal lengths. 
Therefore, when the signals from the two perpendicular arms are combined, the laser phase 
noise in the differenced signals cancels exactly. In space, a laser interferometer gravitational 
wave detector such as LISA [Ij] will have free-flying spacecraft as the end masses, and precise 
equality of the arms is not possible. Other methods must then be used to eliminate laser 
phase noise from the system O, H]. These methods involve a heterodyne measurement 
for each separate arm of the interferometer and data processing that combines data from 
both arms to generate a signal that is free of laser phase noise. In a previous paper ([|J, 
hereafter called paper I), the sensitivity curves for space detectors using these techniques 
were generated by explicitly calculating transfer functions for signal and noise, as modified 
by the data processing algorithms. While the algorithms have been shown ||, in principle, to 
eliminate the laser phase noise in the detectors regardless of the lengths of the two arms, the 
transfer functions have previously only been calculated for the case of equal arms f|, ||, |], [7| . 
In this paper we extend the calculation of the noise and signal transfer functions to the case 
of arbitrarily chosen armlengths. 

One of the goals of paper I was to provide a uniform system for evaluating the sensitivity 
of various configurations of space gravitational detectors. This paper extends that capability 
to configurations in which the armlengths are significantly different from each other. For 
example, a proposal by Bernard Schutz at the 2000 LISA Symposium in Golm, Germany 
f|, suggested a modification to the current LISA design in which a fourth spacecraft is 
inserted in the middle of one of the legs of the interferometer to produce two independent 
interferometers, each with one leg half the length of the other (see Fig. p. The goal of such 
a design was to be able to cross-correlate the independent interferometers to search for the 
stochastic cosmic gravitational wave background. Using the analysis presented here, one 
will be able to determine the sensitivity of such an interferometer and judge the scientific 
value of the proposed modification. 

As in paper I, the analysis begins with the response of a round-trip electromagnetic 
tracking signal to the passage of a gravitational wave, as derived by Estabrook and Wahlquist 
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||. A gravitational wave of amplitude h(t) will produce a Doppler shift Au in the received 
frequency, relative to the outgoing signal with fundamental frequency vo{t). The shift is 
given by 

MW) i 9/ 

= - cos 2w 

v 2 

x [(1 - cos 6)h{t) + 2 cos 9 h(t - r - r cos 9) - (1 + cos 6)h(t - 2r)] , (1) 

where r is the one-way light travel time between spacecraft, 9 is the angle between the line 
connecting the spacecraft and the line of sight to the source, and if) is a principal polarization 
angle of the quadrupole gravitational wave. It is desirable to work in frequency space, so 
h(t) is written in terms of its Fourier transform h(uj). If the Doppler record is sampled for 
a time T then h(t) is related to its Fourier transform by 

\PT r+°° ~ 

h(t) = ^— / h{u)e twt duo , (2) 

where the VT normalization factor is used to keep the power spectrum roughly independent 
of time. Using this definition of the Fourier transform, the frequency shift of Eq. ([!]) can be 
written as 

Au(t, 9, if;) = ^— / - cos(2^) h(u, 9, 0, if;) 



X 



1 - (i) + 2fie- luJT{1 ^ l) - (1 + fi)e~ i2ujT e tujt du , (3) 



where \i = cos^. The quantity that is actually read out in a laser interferometer tracking 
system is phase, so Eq. (||) is integrated to find the phase in cycles 

A(f)(t,6,if;) = J Au(t,9,ip) dt. (4) 

In paper I, a strain-like variable z was formed by dividing the A<f> in Eq. by z/ r and 
the analysis was done using this variable. Since both arms had roughly the same length in 
paper I and carried nearly the same frequency, there was only a scale difference between 
using Acf) and using z as the observable, and linear combinations of z were the same as 
linear combinations of A<f>. However, when the two armlengths are different, this is no 
longer the case, and one must be careful as to what is taken to be the observable for use in 
noise-cancelling data analysis. 

In the laser phase-noise-cancellation algorithms that will be presented in Section II, it 
is relative phase and not strain that can be combined to create laser-noise-free signals. To 



understand how this arises, consider a case where laser signals in two arms are phase-locked 
to each other, with v\ as the frequency of the master laser in the first arm and v 2 = x u i as 
the frequency in the second arm, with x as the ratio of the two frequencies. Then a phase 
noise excursion <50i in the first arm will produce a phase noise excursion x^'Pi i n the second 
arm. Thus it will be linear combinations of Zi = <pij Vi that will allow the two noise terms to 
cancel. Therefore, in this paper, the gravitational wave observable in the i th arm is defined 
to be 



Zi(t,9,tjj) 



T r+°° r n 1 
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dui cos lip h(u) (1-//) + 2ne- WT{1+ »> - (1 + / u)e^ r -e lu)t , (5) 



where Eq. (|3]) has been used to expand Au(t, 9, ip) and where arbitrary constant phases have 
been set to zero in the integration. It should be noted that z% is a different observable than 
the strain variable that was labelled Z{ in paper I. It should also be noted that Zi, as it is 
now defined, has units of time, so Eq. (|5|) gives the time delay in seconds produced by the 
passage of a gravitational wave through the detector. 



II. SENSITIVITY CURVES 
A. Instrument Signal 

Tinto and Armstrong [0 originally showed that the preferred signal for purposes of data 
analysis is not the traditional Michelson combination (difference of both arms), but rather 
a new combination X(t), given in the time domain by fT0| 

X(t) = si(t) - s 2 (t) - si(i - 2r 2 ) + s 2 (t - 2n) 
= Zl (t) - z 2 (t) - Zl (t - 2r 2 ) + z 2 (t - 2n) 

+ni (t) - rn{t - 2r 2 ) - n 2 {t) + n 2 {t - 2t x ) , (6) 

where Si(t) is the data stream from the i th interferometer arm, composed of the signal Zi(t) 
of interest (given by Eq. (^)) and the combined noise spectra in each of the interferometer 
arms, rii(t). The armlengths are taken to be unequal, with armlength Tj in the i th arm. This 
combination is devoid of laser phase noise for all values of the two armlengths T\ and r 2 
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To determine the sensitivity using the X(t) variable, it is necessary to establish a rela- 
tionship between the amplitude of a gravitational wave incident on the detector and the 
size of the X(t) signal put out by the instrument. The noise in the detector will limit this 
sensitivity, and must also be included in the analysis. The part of X(t) containing the 



gravitational wave signal is [ITT 



A(t) = Zl (t) - z 2 (t) - zi(t - 2r 2 ) + z 2 (t - 2n) , (7) 

The transfer function lZ(u), which connects the spectral density of the instrument output, 
Sfr(uj) with the spectral density Sh(uj) in frequency space, is defined via 

S K (u)=S h (u)K(u) , (8) 

where the bar over the A in Eq. (H) indicates an average over source polarization and 
direction. The gravitational wave amplitude spectral density Sh(uj) is defined by 

S h (u) = \h(u;)\ 2 , (9) 

where h(u) is the Fourier amplitude defined in Eq. (fj), so that the mean-square gravitational 
wave strain is given by 

Similarly, the instrumental response S\(u) is defined such that 
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(A^) = — / S K {u)du , (11) 

Z7T JO 

where the brackets indicate a time average. In the next section, the transfer function from 
the gravitational wave amplitude h to the instrument signal A is worked out. 

B. Gravitational Wave Transfer Function 

Let us take the ratio of the two armlengths in the interferometer to be an adjustable 
parameter, /3, taking on values between and 1, such that T\ = t and r 2 = (3r. The average 
power in the part of X(t) which contains the gravitational wave signal is given by 

(A 2 )=lim-/ \A\ 2 dt, (12) 



where A is defined by Eq. ([?]). Using the definition of z from Eq. ([5]) this can be expanded 
to yield 



1 r°° ~ I 
(A 2 ) = — / du h\u)- [TiM + T 2 {u) - 2T 3 (u)} , 



2tt 



(13) 



where 



TAu) 



To(u) 



cos 2 (2ipi) ■ 4 sin 2 (/3u) (l + cos 2 (u) — 2 cos(u) cos(u^i 
2// 1 sin(u) sin{uni) + sin 2 (w 



(14) 



r 3 («) 

with It = CUT, /!. 



(15) 
(16) 



(17) 



cos 2 (2?/'2) • 4sin 2 (w) /i 2 (l + cos 2 ( / 9m) — 2cos(/3-u) cos(/3w//2) 

2/i 2 sin(/3w) sin((3un2) + sin 2 (f3u) , 

cos(2-^i) cos(2 , 02) • 4 sin(u) sin(/3w) r/(w) , 

cos 0j, and where 

rj(u, 9\, 62) = [cos(w) — cos(u^i)] [cos(/3w) — cos(/3u/i 2 )] //1//2 

+ [sin(w) — Hi sin(w/ii)] [sin(/3it) — /x 2 sin(/3u/x 2 )] 

has been defined for convenience. The propagation angles 0j and principal polarization 
angles ipi are defined with respect to the i th arm using the geometric conventions of paper I. 
The expression for the power in the detector, as given by Eq. (|T3|), is a complicated function 
of frequency and of the orientation between the propagation vector of the gravitational wave 
and the interferometer, and represents the antenna pattern for the detector. 

It is customary to describe the average sensitivity of the instrument by considering the 
isotropic power, obtained by averaging the antenna pattern over all propagation vectors and 
all polarizations ]T2]. Using the definition of TZ(u) from Eq. (|S|), with the average isotropic 
power computed using the geometric averaging procedure of paper I with Eqs. (pj] - |i~7[) , the 
gravitational wave transfer function is found to be 



K(u) 



2sin 2 ((3u) 



+ 2shr(u) 



1 + cos 2 (fin] 



cos [u / 

'1 _ 



u- 



■ \ + sin 2 (u) H — - sin(w) cos(m) 



+ sin 2 ((3u) + 



sin(/3w) cos(/9m) 



((3u) 2 J ^ ' ((3u) 3 

— — sin(u) sin(f3u) J de J dfi\ ^1 — 2 sin a) T](u, 81, 6> 2 )j . (18) 

The remaining integral can be evaluated using simple numerical techniques, after relating 

the angular variables as described in Paper I, where: 

sin 7 sine 



sin a 



A. 



(19) 



and 

A*2 = Mi cos 7 + sin 7 cose y 1 — //f . (20) 

Here 7 is the opening angle of the interferometer, and e is the inclination of the gravitational 
wave propagation vector to the plane of the interferometer. The complete gravitational 
wave transfer function is plotted in Fig. ^| for (3 = 1 ("equal arm") and Fig. |] for (3 = 0.01 
( "unequal arm" ) examples. 

As may be seen in the figure, the low-frequency (small u) response of the detector to a 
gravitational wave signal is four orders of magnitude lower for the (3 = 0.01 detector than for 
the equal arm detector, implying that the (amplitude) signal will be two orders of magnitude 
lower - the detected signal level is proportional to the length of the shortest arm. However, 
once the period of the gravitational wave falls inside the light-time of the longest arm, u ~ 1, 
the equal- arm detector ((3 = 1) response begins to fall off while the unequal- arm detector 
((3 = 0.01) response is roughly flat up to a period corresponding to the light-time in the 
shortest arm. 

The dropoff at low frequencies is a result of the fact that the variable X(t) is formed by 
subtracting each Zi from itself, offset by the light-time in the opposite arm. Thus, in the 
low-frequency limit, the two copies of the signal strongly overlap and the signal is almost 
entirely subtracted away. For equal arms, the response of the detector is likewise subtracted 
to zero when an integer number of wavelengths fits in the arm length, as seen in the high- 
frequency portion of the (3=1 curve. For the unequal-arm case, this does not occur, 
because the subtraction of two versions of the signal in each arm are done at different light- 
times in the two arms, so whatever period signal cancels in one arm will not cancel in the 
other. However, as may be seen in the (3 = 0.01 case, the response drops sharply to zero 
at log it ~ 2.5 (equivalent to / ~ 10 5 Hz for LISA armlength of cr = 5 x 10 9 m), where 
exactly one wavelength fits into the short arm and exactly one hundred fit into the long arm. 

However, the response of the detector's X(t) signal is not the whole story. The ability of a 
detector to detect a signal depends on both the signal in the detector and on the competing 
noise. As we shall see in the next section, when the X(t) variable is formed, the noise in 
each arm is likewise subtracted away in most of the places where the signal is lost (e.g., at 
low frequency), so the ratio of signal to noise remains high. 
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C. Noise transfer function using the X(t) variable 



The noise sources for LISA may be divided into categories in two different ways. First, 
a noise source may be either one-way (affecting only the incoming or the outgoing signal 
at a spacecraft, but not both) or two-way, (affecting both incoming and outgoing signals at 
the same time). A one-way noise source will have a transfer function of 2, since there are 2 



spacecraft in each leg contributing equal amounts of such noise |T3|]. The transfer function 
for two-way noise sources, however, will be more complicated due to the internal correlation. 
A single two-way noise fluctuation in the central spacecraft of the interferometer will affect 
the incoming signal immediately, and then, a round-trip light-time later, will affect the 
measured signal again in the same way. In the time domain, the effect in the i th arm of a 
fluctuation n(t) will be rii(t) = n{t) + n(t — 2rj). The transfer function for this time-delayed 
sum is 4 cos 2 (2tt f Ti) . If an end spacecraft has noise that affects both incoming and outgoing 
beams, it will affect them at almost the same time, with no delay, giving a transfer function 
contribution of 4. The noise transfer function for a single arm for a two-way noise source is 
therefore 

4 + 4cos 2 (2vr/r i ) . (21) 

Examples of one-way noise are thermal noise in the laser receiver electronics or a mechanical 
change in the optical pathlength in the outgoing laser signal before it gets to the main 
telescope optics. Examples of two-way noise are parasitic forces on the accelerometer proof 
mass or thermal changes in the optical pathlength in the main telescope. 

A second way in which noise sources may be classified is by how they scale when there 
is a change in armlength in the interferometer. The first type of noise in this classification 
scheme is what we call "position noise", in which the size of the noise in radians of phase 
is independent of the length of the arm. Accelerometer noise and thermal noise in the 
laser electronics are examples of position noise. The second type of noise is what we call 
"strain noise", in which the size of the noise scales with armlength. Examples of strain 
noise include shot noise and pointing jitter (if it is dominated by low power in the incoming 
beacon). Position noises may be either one-way or two-way, but we can think of no two-way 
strain noise sources. 

The transfer functions that connect the noise in the instrument to the X variable depend 
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on the type of noise. We begin by considering the noise terms in Eq. (P): 

a(t) = m(t) - n 2 (t) - [m(t - 2(3t) - n 2 (t - 2r)] . (22) 

We then go to the frequency domain, squaring and time-averaging to obtain the power 
spectrum. 

(a 2 ) = — f du 4 \h\ sin 2 ((3u) + h\ sin 2 (w)l , (23) 

2,71 J L - 1 

where cross-terms (e.g., rli^) have been neglected under the assumption that noise in the 
two arms will be independent and uncorrelated. Note that n\ is the power spectrum in the 
long arm (length r) and n\ is the power spectrum in the short arm (length f3r). 

Since the noise in the detectors includes different types, with different transfer functions, 
it is not possible to write a single transfer function giving the response of the X variable 
to noise, so let us consider the various noise categories one at a time. We first consider 
position noise, for which n 2 = n\ = n 2 . Then, using Eq. (E^), we find the transfer function 
for one-way position noise to be 

^i = 8 (sin 2 (/?w) + sin 2 (w)) , (24) 

where, as we noted above, there is a factor of 2 representing the noise from the two spacecraft 
in each arm. Two-way position noise must include the transfer function from Eq. (pi]), giving 



Uo = 16 



sin 2 (/3w) (l + cos 2 (w)) + sin 2 (w) (l + cos 2 (pu))] . (25) 



Strain noise scales with armlength, and is hence smaller in the shorter arm, so that h 2 = 
h\ = n\/ (3 2 . Its transfer function is therefore 

K s = 8 (sm 2 (pu) + p 2 sm 2 (u)) , (26) 

where the factor of 2 for the two spacecraft has again been included. 

When (3=1, the transfer functions for strain noise and one-way position noise (Eqs. |4] 



and |26|) are identical and have zeros at u n = nir, where n is zero or a positive integer. These 
are exactly the places where the (3 = 1 transfer function for gravitational wave signal (Fig. 
H) has its zeros. When (3 < 1, the situation is more complicated. Both 7Z± and 1Z S share the 
sin 2 ((3u) term which will go to zero at u = and at multiples of u = ft/ (3. The sin 2 (it) terms 
in TZi and 1Z S have their zeros at multiples of the lower frequency, u = ix. In 1Z\, this term 
will be larger than sin 2 ((3u) term at low frequencies, since near u = 0, sin 2 (u) ~ u 2 , while 



sin 2 (/?-u) ~ j3 2 u 2 . In 7Z S , these terms will be equal in the low-frequency limit, because of the 
factor f3 2 that multiplies the sin 2 (w) term. Thus, in the low frequency limit, the strain noise 
transfer function will be 2f3 2 times the the one-way position noise transfer function. When 
/3 < 1, the transfer function for one-way position noise will have sharp drops at multiples 
of u = 7r, down to the level of its sin 2 (flu) term. These behaviors are shown in Fig. |] and 
Fig. § 



D. Sensitivity curve 

The signal to noise ratio is the ratio of the signal power in the detector to the noise power 
in the detector: 

SNR = s s n s + + s 2 n 2 ■ (27) 

where S s , Si, and S 2 are the spectra of strain noise and one-way and two-way position noise, 
respectively, and 1Z is the gravitational wave transfer function given by Eq. fllBp . Setting 
SNR = 1 and solving for hf = y/Sh yields the instrument sensitivity curve as defined in 
paper I: 

, /7T \s 8 n 8 + S&! + S 2 K 2 

hf = V S h = )J ^ • (28) 

where TZ is the gravitational wave transfer function, given by Eq. (|18|). 

Figures ^ and [7] show the sensitivity curves, computed using Eq. (P%D, for /3 — 1 and 
(3 = 0.01 respectively. The noise values used are taken to be the LISA target design values 
(computed as described in paper I). The shot noise and acceleration noise levels are set at 
the standard LISA values. In addition, a flat one-way position noise spectrum is assumed 
at l/10 th the LISA shot-noise value. Also plotted in Figures^ and are sensitivity curves 
representing each of the three components of the total noise, taken one at a time. 



III. DISCUSSION 

As may be seen in Fig. |7|, the low-frequency sensitivity for unequal arms, being set by the 
two-way position noise in the accelerometer, is degraded over the equal-arm case by the ratio 
of the two arms. In other words, the sensitivity at lowest frequencies is set by the sensitivity 
of the shortest arm. At middle and high frequencies, the situation is more complicated. If 
the dominant noise is strain noise, then the sensitivity is independent of j3 in this frequency 
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range. However, if the dominant noise is position noise, then the sensitivity curve at high 
frequencies will rise in proportion to /3, though its flat floor will extend to higher frequency, 
from the 1/(2ttt) of the equal-arm case to 1/(27t/3t) when the armlength ratio is /3. 

The implications of these results for mission design are obvious. If the armlengths are not 
equal, the low-frequency sensitivity is degraded by a factor 1//3, the ratio of the armlengths. 
If the high-frequency noise can be guaranteed to be strain noise, even in the shorter arm, 
then the high-frequency sensitivity is unaffected by the unequal arms. If the noise at high 
frequency is dominated by position noise, then the high frequency sensitivity is degraded by 
the factor 1/(3, but the sensitivity remains flat up to a frequency l/(27r/3r), where it turns 
over and joins the strain noise curve. Thus, as long as the position-noise sources can be kept 
well below the shot noise and other strain-noise contributions, a change in armlength ratio 
from strict equality will not degrade the high-frequency portion of the sensitivity curves. 
However, as the length of one of the arms is shortened, small position noise sources will 
become important and eventually dominate. 

Let us consider the example of Schutz's 4-spacecraft configuration (Fig. |l|). Since this 
configuration will have (3 = 0.5, the low-frequency sensitivity curve will be a factor of 2 
higher (hence less sensitive). The current error budget for LISA assumes that the high- 
frequency portion of the window is dominated by position noise approximately three times 
the shot noise. If this remains the case, then the high-frequency section of the curve will 
likewise be a factor of 2 higher up to a frequency twice as high as the LISA sensitivity "knee" 
at / = 1/(27tt), at which point it would turn up and join the current LISA high-frequency 
ramp. The shot noise is determined by the power of the laser and by the size and efficiency of 
the optics, and there is nothing beyond brute-force improvements in these parameters that 
will lower the shot noise. The contributions to position noise, on the other hand, are due to 
optics quality, the attitude control system, Brownian noise in the electronics, thermal noise 
in the optical path length, etc. These are more complex and are amenable to reduction by 
careful or innovative engineering design. If these noise sources can be reduced to a fraction of 
the shot noise, not only will the LISA noise floor be reduced by a factor of 4, but the Schutz 
modification will have high-frequency performance that is undiminished by the reduction of 
the length of one arm. 

Finally, we describe a totally unfeasible mission design that is nevertheless interesting 
for instructive purposes. Let us consider a two-spacecraft "interferometer", where one of 
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the spacecraft contains a fiber optic delay line, of length 5 km, that acts as the second arm 
of the interferometer. If the distance between the two spacecraft is 5 x 10 6 km, we will 
have (3 = 1CT 6 . The use of the X(t) variable will eliminate laser phase noise, exactly as it 
does in arms that are more nearly equal. A rigidly-attached reflector at the far end of the 
fiber-optic line would eliminate accelerometer noise, but, of course, would replace it with 
thermal fluctuation in the optical path length in the fiber. However, a concatenation of 
fibers with well-chosen thermal pathlength coefficients could produce a fiber tuned to have a 
coefficient very near zero. This, combined with multilevel thermal isolation, could keep this 
noise source very small. The key to the sensitivity of this configuration is the position noise. 
If a way can be found to reduce position noise to less than 10~ 6 of the LISA shot noise, 
then this two-spacecraft interferometer would have the same sensitivity as a conventional 
three-spacecraft interferometer. 
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FIG. 1: An unequal arm geometry used here assumes two arms of length r and (3t, with an enclosed 
angle 7 (the interferometer opening angle). Depicted here is the nominal LISA constellation of three 
spacecraft in an equilateral triangle, and a proposed extension which places a fourth spacecraft 
midway down one of the arms. 

FIG. 2: The dimensionless gravitational wave transfer function, TZ/t 2 , plotted against the dimen- 
sionless frequency parameter u = ojt, for value of (3 = 1.0. 

FIG. 4: The noise transfer functions for (3 = 1 as functions of the dimensionless frequency 
parameter u = ojt. Notice that the transfer function for position noise (TZi) is identical to the 
transfer function for shot noise (7£ s )in the (3 = 1 limit. 



FIG. 5: The noise transfer functions for (3 = 0.01 as functions of the dimensionless frequency 
parameter u = lot. 



FIG. 6: The sensitivity curve (SNR = 1) for (3=1. Overlayed are the sensitivity curves for each 
of the individual noise spectra (acceleration noise, shot noise, position noise). The noise spectra 
are taken to be at the LISA target design values, except position noise, which is taken to be l/10* h 
the LISA value. 



FIG. 7: The sensitivity curve (SNR = 1) for [3 = 0.01. Overlayed are the sensitivity curves for 
each of the individual noise spectra, as in the previous figure. 

FIG. 3: The dimensionless gravitational wave transfer function, TZ/t 2 plotted against the dimen- 
sionless parameter u = ojt, for value of (3 = 0.01. 
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